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Abstract 

In the standard setting of one-way AN OVA with normal errors, a new algo- 
rithm, called the Step Down Maximum Mean Selection Algorithm (SDMMSA), 
is proposed to estimate the treatment means under an assumption that the 
treatment mean is nondecreasing in the factor level. We prove that i) the SD- 
MMSA and the Pooled Adjacent Violator Algorithm (PAVA), a widely used 
algorithm in many problems, generate the same estimators for normal means, 
ii) the estimators are the mle's, and iii) the distribution of each of the esti- 
mators is stochastically nondecreasing in each of the treatment means. As an 
application of this stochastic ordering, a sequence of null hypotheses to iden- 
tify the minimum effective dose (MED) is formulated under the assumption of 
monotone treatment (dose) means. A step-up testing procedure, which controls 



the experimentwise error rate in the strong sense, is constructed. When the 
MED=1, the proposed test is uniformly more powerful than Hsu and Berger's 
(1999). 



Some key words: Closed test method; Experimentwise error rate; Maximum 
likelihood estimator; Step-up tests. 

1 Introduction. 

A situation frequently encountered in dose-response studies is identifying the mini- 
mum effective dose [MED). The MED is defined as the lowest dose such that the 
mean response is better than that of a zero-dose control by a clinically significant 
difference. Finding the MED is important since high doses often turn out to have 
undesirable side effects. 

Consider the one-way layout model 



for i = 0, ■■■,k,j — 1, ...,ni, where //j's are the unknown response means at different 
dose levels and ~ ^'"(O, a^) are the independent errors with an unknown variance. 
The parameter space is 



(here, for simplicity, we omit cr in if), and the sufficient statistics are the sample 
means, Yi, and the mean squared error, denoted by S'^. Assume i = is the control 
group. One goal is to find the smallest positive integer satisfying /in > /iq + S for a 
chnically significant difference constant 6 >0. We call A^ the minimum effective dose 
(MED). Determination of the MED usually is done by step-down test procedures, 
see Williams (1971), Ruberg (1989), Tamhane, Hochbcrg, and Dunnett (1996), Hsu 
and Berger (1999), and Hellmich and Lehmacher (2005), among others. Tamhane, 




(1.1) 
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Hochberg and Dunnctt (1996) indeed proposed a step-up procedure SUIP to iden- 
tify the MED. The SUIP procedure is based on the step-up procedure of Dunnett 
and Tamhane (1992), which controls the experimentwise error rate only for balanced 
designs. However, Dunnett and Tamhane (1995) 's step-up procedure for unbalanced 
designs case cannot control the experimentwise error rate. Liu (1997) proposed a 
method of calculating the critical values of the step-up procedure by Dunnett and 
Tamhane (1995). The SUIP procedure does not make use of the monotonicity, there- 
fore its power should not be high. Intuitively it seems that step-down procedures infer 
a larger dose as the MED. Therefore, it is of interest to have a step-up procedure to 
use the monotonicity to increase its power. 

To derive a test of level-a, one needs to find an appropriate statistic and its least 
favorable distribution in the null hypothesis. Thus a stochastic ordering for the test 
statistic is needed. The desired statistic, the estimator of //j, should be: a) non- 
decreasing in i, b) and is also nondecreasing in each of Y/s. The PAVA algorithm 
generates the estimators that achieve a). However, it is difficult to show b) directly 
for these estimators using the PAVA. The PAVA was first proposed by Ayel, Brunk, 
Ewing, Reid and Silverman (1955), and was introduced to estimate the monotone 
proportions in independent binomial experiments. Surprisingly, it has many appli- 
cations in normal, Poisson and multinomial distributions, etc. See more details in 
Robertson, Wright and Dykstra (1988). The PAVA is an iterative algorithm, each 
step is very simple to implement, however, it does not have a closed form for the final 
estimator. Hence, it is hard to establish analytic properties for the estimator. Notice 
these, a new algorithm, the SDMMSA, is proposed to overcome the drawbacks. We 
will show that the two algorithms yield the same estimators and each estimator is 
a monotone function of each The second fact is critical to determine the least 
favorable distribution in the null hypothesis space. 
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The rest of the article is organized as follows. Section 2 provides a new algorithm 
to generate estimators for /ij's and discusses their analytic properties. In particular, 
a stochastic ordering for the distributions of the proposed estimators is established. 
In Section 3, one application of the stochastic ordering established in Section 2 is 
given to identify the MED. A step-up multiple test procedure that controls the ex- 
perimentwise error rate in the strong sense is provided by constructing a sequence 
of increasing rejection regions of level-a for each null hypothesis in f l3.25p . and the 
proposed procedure is illustrated on a real data set. Section 4 concludes with some 
discussion. 

2 A new algorithm to construct the estimator of 
fji and some analytic results. 

In this section, an estimator of /i,, denoted by fii, for any integer i E [1, fc] under H 
is first constructed iteratively. Then three facts are established: fii is the same as the 
estimator generated by the PAVA; fii is the mle under H, and the distribution of fii 
is stochastically non-decreasing in each fij. 

2.1 A new iterative algorithm to construct fii. 

Let 

n,,, = n,, Yi, = yi<t<j<k (2.3) 

h=i ^hj 

be the sample size and the sample mean of a combined sample of treatments i through 
j, respectively. 

Step 1). We construct fii starting from i = k using the data set {(Yi, ni)}f^^. Let 
Ai = {j : Yj^k = max{i<j'<k}{Yj'^k}} (2.4) 
be a subset of {1, k}{Ai contains a single element with probabihty one), and let 

ii = min{Ai}. (2.5) 
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Then 

/i/=^F,,,,, V «G (2.6) 

If ii = 1, then all /ij's are defined and stop; otherwise go to the next step. Step 2). 
Note in this step ii — 1 < k — 1. Repeat Step 1 but using the data set {(Fj, n,)}*!"^"^. 
i.e., let 

^2 = {j ■■ yj,ii-i = mO'X{i<j'<ii-i}{^j',ii-i}} (2-7) 
be a subset of {1, — 1}, and let 

^2 = min{A2}. (2.8) 

Then 

ixi=^Y,,^i,^i, V tE[t2,ii-l]. (2.9) 

If ^2 = 1, then all /xj's are defined and stop; otherwise repeat this process for a 
number of times, say h times, until ih = 1- Such an integer h exists, because ij 
strictly decreases in j. Since i/^ = 1, then all defined and the construction 

on /ij's is complete. We name this the step-down-maximum-mean-selection algorithm 
(SDMMSA). 

Remark 1. There exists partition, u'^^i[iu, iu-i — 1], for {1, k} with ~ 1 k. 
Following the construction of fii, each fii is the sample mean of a combined sample 
of treatment (s) belonging to the interval in the partition that includes treatment i. 
Also fii is constant in i on each integer interval [iujiu-i ~ 1] for u = 1, ...h, as shown 
in (12. 6p and (12.91) . and /tj is strictly increasing when i moves from [iujK-i — 1] to 
[iu','iu'-i — 1] for u > u', as shown in (12.41) and (12.71) . Therefore, fii is nondecreasing 
in i for z G [1 , A;] . □ 

Lemma 1 For partition ^u=i[^u,'i-u-i ~ 1] given in Remark 1, Yi^_i < Yi^ for any 
u e [1, h]. 
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Proof. Since the SDMMSA repeats itself in each step, without loss of generality, we 
only need to prove Lemma [1] for u = 1. i.e., Fj^-i < Vij. 

Suppose this is not true, i.e., Fj^ < Yi^^i. Note Yi^^k > '^h+i,k by the definition of 
ii, then 

Fn > Fi+l,fc- (2.10) 

Similarly, < Yi^^k is true due to Yi^-i^k < Fi.fc- Therefore, Yi^ < Yi^-i < Yi^^k, 

which implies 

Fi ^ Fi+i,fc) 

a contradiction to ( I2.10p . □ 

Remark 2. If Fj is nondecreasing in i G [1, k], then the partition, U'^^-^liu, iu-i — 1], 
for {1, k} given in Remark 1 satisfies i) F is constant when i G [iu,'iu-i — 1], and 
ii) is strictly increasing when i moves from [iujiu-i — 1] to [iu',iu'-i — 1] for u > u'. 
Therefore, /ij = Fj for z G [1, /c]. □ 

Example 1. Consider the data in Table 1, taken from Ruberg (1995). There are nine 
{k = 9) active dose groups and a zero dose control group with six (n^ = 6, i = 0, 9) 
animals/group in the experiment. Following Step 1, we obtain ii = 9 and then 
fig = Fg; following Step 2, we obtain ^2 = 6, then fi^ = fi^ = fiQ and is equal to 
F6,8 = 73.77, the sample mean of the combined sample for i = 8, 7, 6. The construction 
of all /ij's ends at Step 7(=h) and their values are reported in Table 1. The partition 
given in Remark 1 is now 

[9]U[6,8]U[5]U[4]U[3]U[2]U[1], 

with a notation of [i] = □ 
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2.2 The relationship between /i^, /i™'^ and /if""'*. 

So far, the estimator of fii under H typically is obtained following the pooled- adjacent- 
violators algorithm(PAVA, described later), for example, see Barlow, Bartholomew, 
Bremner and Brunk (1972), Robertson, Wright, and Dykstra (1988), and SilvapuUe 
and Sen (2005). Denote this estimator by fi^"'^"'. Now we show that /t™'^ = /^i = /tf^^" 
in Theorem [1] and Theorem [2] below. 

Theorem 1 Let be the maximum likelihood estimator for fii under H . Then 
fii = fif^, Wie[l,k]. Therefore, fn = fif^, Vi G [1, fc] . 

Proof of Theorem [H Taking the log transformation on the joint pdf of Yij, it is 
easy to see that /i™''^ minimizes 

k h h 

1=1 u=l j=iu u=l 

under H, where the intervals iu-i — 1] for u = 1, h are given in Remark 1. 

Now focus on each /j. Without loss of generality, focus on /i, then 

k _ k _ 

/i(/ii,...,/Xfc) = nj{Yj~iijf = n,[{Yj-fijY + 2{Yj-fLj){fLj-fij) + {fij-fijy]. 

j=il j=h 

Rearrange the terms above and note fij = Yi-^^k for j G [ii, k], then 

k _ k _ 

/i(/ii, ...,/Xfc) = {Y nj[{Yj - fijf + {fij - lijf]} + 2 ^ rijiXj - =^ h + h- 

j=il j=ii 

It is obvious that Ii is minimized at fij = fij for j G [ii, k]; for I2, apply Abel's partial 
summation formula and obtain 

ii il+l 

/2 = 2 ^ rijifij - Yj)fij = 2 ^ djifij - fij^i) + di^fii^ = 2 ^ djifij - fij_i), 

j=k j=k j=k 

where dj = Y^i^kn^ifiv — Y^) > and di-^ = due to the definition of ii. Also note 
fij > ^j-i. Thus I2 is nonnegative and achieves its minimum at /ij^ = ... = fik- 



Therefore, combining Ji and J2, we conclude /ifc) is minimized at fij = fij 

for j G [ii,k]. 

Repeat the same argument on /2 through fh, each fu is minimized at fij = fij 
for j G [iu^iu-i — !]• Lastly, note fij nondecreasing, we conclude fi^^'^ = fij for any 
J e [l,k]. □ 

Theorem 2 For any z G [1, k], 

/2.(Fi,...,n) = /ir'^(?i,...,n). (2.11) 

For a data set of '^i)}i=i) ^^e PAVA proceeds as follows: 
Step 0-PAVA). If Yi is nondecreasing in i for i G [1, k], then fi^"'^"' = Yf, otherwise, 
go to the next step. 

Step 1-PAVA). Pick any consecutive pair {Yj,Yj^i) with Yj > y,+i, let ji be the 
smallest integer so that Fj = Yj for i G [ji,j] and let j„ be the largest integer so that 

Yi = Yj+i for i e[j + l,ju]- For each i G [ji,ju], replace each Yi by l^j,,j„(= ^^jj? ), 

and then obtain a new data set of {{ai,ni)}i^^, where Oj = Yi for i ^ [j^ju] and 
tti = Yj^j^ for i G [ji, Jm]. Note two facts: Yi is non-increasing for i G [j;, ju], and the 
number of different Oj's is strictly less than that of Fj's. 

Step 2-PAVA) Repeat this process on {(aj, ni)}i^i until all aj's are nondecreasing. 
Then fi^"'^"' = ai. Since the number of different a^'s is strictly less than that in the 
previous step, this algorithm has to stop in a finite steps. □ 

Proof of Theorem [2]. When Yi is nondecreasing in i G [l,k], then (12. lip is true 
due to Step 0-PAVA) and Remark 2. When Yj > Y^+i for some j, let and [ji,ju] 
be given in Step 1-PAVA). It suffices to show 

f,i{Yi, Ffc) = P^i{ai, flfc), (2.12) 
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for any i G [l,k]. Let U'^^iliu, iu-i — 1] be the partition of [1, A;] given in Remark 
1 using data {(Yi,ni)}^^^. The integer j„ has to belong to one of these intervals in 
the partition, say [iuj,'iuj — !]• Since Yi is non-increasing on [ji,ju] as shown in Step 
1-PAVA), by Lemma [H i^^ < ji. Thus [ji,ju] is a subset of [iup^uj ~ l]? an interval 
in the partition. Let ^u=i[^'u^^'u-i ~ 1] be the partition of [1,A;] given in Remark 1 
but using data {{ai.rii)}^^^. Therefore, [ihiv\ also has to be a subset of one of these 
intervals. 

Case I). If [ii,iv\ C [«i,A;], i.e., Uj = 1, consider 

def J2u=i ''^uQ'u 

for i G [1, k]. Note Oj^fc = Yi^k for i ^ [ji + l,ju] and i[ < ji, then i[ = ii. Therefore, 

fXi{Yi, Yk) = = — ^ = /i^ai, Ofc), 

for any i G [ii,/:;]. For any i < ii, since fii only depends on Yi through Fj^_i(or ai 
through Oj^-i) and Yi = ai, we conclude (12.120 . 

Case II). If [ji,ju] C [i2,ii — 1], i.e., Uj = 2, we only need to show i[ = ii. Then, 
similar to Case I) above, (12.121) is established. To prove i[ = ii, first note 

for any u > ii. So 

i[ < ii (2.13) 

by the definition of i'l- 

Suppose i'^ < ii, we will construct a contradiction. Note i'^ ^ Uhju\ because Oj is 
non- increasing on (j;,j„]. Thus, 

a contradiction with the definition of ii. Hence i'^ > ii. Together with (I2.13p . we 
conclude i[ = ii. 

9 



For the other cases of uj = 3, h, similar to Case II), we can show i„ = i'^ for all 
u = 2, h. Hence, h' = h and two partitions, U'^^-^[iu, iu-i — 1] and U'^=i[iu, i'u-i ~ 1] 
are identical. Also note that [ji,ju] is contained in one interval [iuj,'iuj-i — 1], (12.121) 
is established. □ 

Remark 3. Although fii and fi^"'^"' generated by two algorithms are identical, there 
are several advantages of the SDMMSA over the PAVA. First, it is clear from the 
SDMMSA that fii is uniquely defined, but not clear for /tf"^" from the PAVA, since the 
latter needs to show /tf"^" must be the same no matter where to start the algorithm, 
which is not obvious at all. Secondly, fij has a closed form, if j G — 

1] for some jo, where ij^ is given in Remark 1, while /if"^'^ does not. This fact is 
important for establishing the stochastic ordering of fii as shown in the next section. 
Thirdly, it was mentioned, for example, in Robertson, Wright and Dykstra (1988, 
p. 10) that jl^"'^"' is the mle. To the best knowledge of the authors, no rigorous proof 
has been given. With fii, we proved fii = fi^^^ and fii = fif"^"" . Thus fFj^^"" = fi^'"'^. 
Lastly, regarding the computation, the SDMMSA is easier to code than the PAVA 
since at each step of the SDMMSA a certain number of the final estimators {fii) are 
defined. □ 

2.3 A stochastic ordering of fii. 

We provide another major result in this paper which establishes a stochastic ordering 
for each fii in terms of each of /ij's. Let 

fl^,k = f^^{Yl,...,Yk) (2.14) 

be the estimator of /ij obtained from the sample {(Fj,nj)}jL^ following Steps 1 and 2 
in Section 2.1. So the distribution of fii^k depends on /ii through /Xfc and a. 
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Theorem 3 For each i and j in [1,A;], fii^k, a function ofYj, is nondecreasing 
when the other Yji 's are held fixed. Therefore, fii k is stochastically nondecreasing in 
fij when the other ^ji's are held fixed, i.e., P{fii^k > x) is a nondecreasing function 
of 11 j for any real number x. 

Proof of Theorem [31 We will prove the monotonicity of fii^k in each Yj by induction 
on k. 

First for the case of A; = 1, fii^i = Yi is nondecreasing in Yi. 

Assume that, for the case of k = m, fii^m is nondecreasing in Yj for any i and j no 
larger than m. Now consider the case of = m + 1. Following Step 1, ii depends on 
k, so write it as ii(fc), i.e., obtain ii{k) using (Fi,ni) through (Ffc,^^). Claim 



m -\- \, if Ym+l > /iji(m),m) 

il(m), if Y^+l G (/ii2(m),m5 Aii{m),m] 



ii{m + 1) 



= < 




(2.15) 



When Yr 



m- 



^ ih{m), if Y^+i < fii^ 
+1 > /in(m),m(= /im,m), for any j E [l,m], note 




then ii{m + 1) 



m + 1. 



When Ym+l G {fii2(m),m, f^him),m] = (^^ia (m),n {m)-l , (m),m] • i) For j G [ii{m) , TU + 

1], note Uj^rn+i - ^ii{m),m+i < 0, then 



> 




,m 
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and conclude ii{m + 1) < ii(m). 
ii) For j G [l,ii{m) — 1], note 

-r> ^j,ii(m) — l'^j,ii{m) — l ~l~ ^ii(m),m^i(m),m ~l~ ^m+l^m+l 

and l^-,ii(m)-i < Yi2(m),ii(m)-i < then 

^l(m),m+l ^',m+l 

'^j,n(m) — l^ii(ni),m^i(m),m ^~ ^j/,ii(m)— l''^m+l •^^+1 ''^j,ii(m,) — 1 ('^ii(m),m ^~ '^m+l)^',ii(m) — 1 

^ii (m),m+l'^j,m+l 

> 0, 

and conclude ii{m + 1) > ii{m) — 1. Therefore, combining i) and ii) we obtain 
ii(m + 1) = ii(m) when Y^+i G {fii^{m),m 

For the other cases of Y^+i, (12.151) can be established in a similar way. There- 
fore, we conclude that fii^m+i depends on {{Yj,nj)}Y=i through {(Aj.m, "^jOljli and 
{Ym+i,n,n+i). So write 

Also write fii^m+i as 

(F,), (2.17) 

since the other Yji are fixed. We will use any one of the above two notations whenever 
it is convenient. For y < y', we need to show the monotonicity below 

P'i,m+liy) < P'i,m+liy'), (2.18) 

which establishes the theorem, in the following two cases. 

Case 1: j = m + 1. Since fij'^m does not involve Ym+i for all j' < m, (I2.18P changes 

to 

(2.19) 
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which is estabhshed in Lemma [2] by noting /tj^m is nondecreasing in i < m. Case 2: 
j < m + 1. Since (If^rniy) < i^j',m{y') for all f < m by the induction assumption on 
the case of k = m, f l2.18p changes to 

which is established in Lemma [31 Therefore, the proof of the monotonicity of fii^k in 
each Yj is complete. 

Since 1^'s are independent random variables, and each is stochastically increasing 
in fij, fii^k, as a nondecreasing function of each Yj, is also stochastically nondecreasing 
in fij. See, for example, Alam and Rizvi (1966) or Lemma 2 in Wu and Wang (2007). 
The proof of Theorem [3] is complete. □ 

Lemma 2 Let fii^rn+i{Yi, ■■■,Ym, y) be the estimator following Steps 1 and 2 on a date 
set {{Y^,n^)}'^^i with Ym+i = y- Then 

{Yu...,Y^,y)< P'i,m+l(Yi, Ym, y ) 

if Yy is nondecreasing in v E [1, m] and y < y' . 

Proof of Lemma [2]. Now write ii{m + 1) introduced in (12.151) as ii{y)- 
First note «i(y) < ii{y'), which follows (I2.15P and y < y'. 
Secondly, claim 

■ / ^ . r r, -,1 -f> ^ -fv T denoted by . 

ti{y) = Tmn{v e[l,m + l]:Y,>Yy^^+i} = ^ A. (2.21) 

Note 

V V — ^'"0^'" ~ ^,m.+ l) _ ^viYy — Yy^l^m+l) /r) r)r)\ 

J-v,m+l ~ J-v+l,m+l — — • [Z.ZZ) 

Therefore, Yy^m+i is nonincreasing in v when v > A. Hence ii{y) < A. It is obvious 
that YA-i^m+i < ?A,m+i following the first equality in (I2.22p . Since Y^ nondecreasing 
in f G [1, m], Y^^m+i is nondecreasing in v when v < A — 1. Thus (I2.2ip is established. 
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Thirdly, a) when i < ii{y), both fii^rn+i{Yi, ...,Ym,y) and ...,Ym,y') are 

constructed based on {Y^y^^i^~^, a subset of {F^jJJL]^ which is nondecreasing in v. 
Thus jli^rn+i{Yi, Ym, y) = % = /ii,m+i(Fi, F„, y') following Remark 2. 

b) When i E [ii{y) , ii{y') - 1], F„, y) = Yi^(y),rn+i with = y; 
while /ii_m+i(Fi5 •••5 ^m, 2/') = > Yi-^{y),m+i following fl2.2ip and Y^ nondecreasing in 
V G [1, m]. 

c) When i G [ii{y'),m + 1], F^, ?/) = with Ym+i = y; 
while fli^rn+iiYi,...,Yrn,y') = F^(y),m+i with = y' . Then /ii,^+i(Fi, F„, < 
fii^rn+i{Yi, ...,Ym,y') due to zi(y) < ii{y'), Y^ nondecreasing in G [l,m], y < y' and 
fl2.2ip . The proof is complete. □ 

Lemma 3 Let fii^m+iiXi, ...,Ym,Ym+i) be the estimator following Steps 1 and 2 on a 
date set {{Y^.n^)}^^^ . Then 

/'■i,m+l(F) •■■) Ymi Fn+l) ^ Ai,m+l(F; Fn? Fn+l) 

z/Yy and Y^ are both nondecreasing and Y^ < Y^ for v G [Ij'W']. 
Proof of Lemma [Si Claim 

/^i,m+l(F_? •••) F)— 1? F)) F)+l") Fti) Fn+l) ^ /'■ijm+l (F_ ; F)— 1? F;' Fn? Fn+l) 

(2.23) 

for any v G [1, m] if F < ... < Y,_i < K < Fj < F^+i < ... < If (E23]) is true, 
then 

To show fl2.23l) . now write ii{m + 1) introduced in fl2.15p as ii(Ft,)(note ii{y) 
introduced in the beginning of the proof of Lemma [2] has a different argument y = 
Ym+i). Following (lOB and the second equality of ^M), «i(F,) < Similar 
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to the proof of Lemma [2], (12.231) can be shown in three cases a) i < ii{y), b) i G 
[ii{y) , ii{y') — 1] and c) i E [ii{y'),m + 1], and the proof is complete. □ 

In short, in this section, we proposed the SDMMSA to generate estimators for 
monotone normal means, fii, proved that the SDMMSA and the PAVA are equiv- 
alent, both generate the mle's, and the distribution of the proposed estimator is 
stochastically nondecreasing when fii goes larger. The last is to be used to derive a 
test to detect the MED in the response-dose study as shown in the next section. 

3 A step- up testing procedure to detect the MED. 

Now return to the problem of finding the minimum effective dose(MED). First, we 
formulate this as a multiple test problem by proposing a sequence of decreasing null 
hypotheses. Then a general result that identifies the least favorable distribution is 
provided. Finally, a sequence of increasing rejection regions of level-a is constructed. 

3.1 Motivation. 

Let 

Xi = Yi- Yq and r]i = fii - /iq for i = 1, k. (3.24) 

Since the MED is to be found, one should start the search from i = 1 instead of 
i = k. Therefore, a step-up procedure seems more reasonable than a step-down one. 
To establish = 1, some authors (see, for example Hsu and Berger (1999)) compare 
min{Xj : j > 1} with 6 and claim A^ = 1 if min{Xj : j > 1} — 5 is large in the 
unit of S. Roughly speaking, they use min{Xj : j > i} to estimate rji. This does 
not fully utilize the assumption of the monotonicity on means. So we propose using 

def „ - 

the maximum likelihood estimator of rji, denoted by fji = fii — Yq, as a test statistic, 
where fii is given by the SDMMSA in Section 2.1. If 171 —5 is larger than a multiple of 
S", then claim A^ = 1 and stop; otherwise compare fi2 — S with S. Repeat this process 
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until we find an N so that tin — S is much larger than S. If no such N can be found, 
then the MED does not exist. 
To identify A^(MED), let 



be the set of null hypotheses of interest in this section. For each i >1, the alternative 
Hai claims r]i > S. If a certain Hai is established, then N <i due to the monotonicity 
in for i > 1. Therefore, N should be equal to the smallest i so that HAi is true. 
For the strong control of the experimentwise error rate, it is clear that i^oi' is a subset 
of if i < i' due to the monotonicity(i.e., H^i is decreasing). Therefore, C itself is 
closed under the operation of intersection. The closed test procedure (Marcus, Peritz 
and Gabriel, 1976) can be applied on C to construct a step-up testing procedure with 
the experimentwise error rate controlled at a in the strong sense (see, for example, 
Hochberg and Tamhane (1987) for a definition) as long as a level-a test is constructed 
for each iJoi- Let Ri be a rejection region for Hoi for any i between 1 and k. In order 
to strongly control the experimentwise error rate, as well as being powerful, region 
Ri should satisfy the following two properties: 
*) Ri is of level a, i.e., sup ^^fj^.Pfj_{Ri) — a. 

**) Ri is increasing in i. i.e. Ri C Ri' if i < i'. Thus Ri — ^\^^Hg^,cHoiR^' — ^i'=iRi'- 
If these two are satisfied, then the multiple tests, which assert Hai if and only if Ri 
occurs, strongly control the experimentwise error rate at level a, which is the main 
result of this section. 

3.2 A general result. 

Theorem 4 Let T{ti, tk) and gi{ti, tk) for i — 1, k be non-decreasing func- 
tion for any ti when the other tj 's are held constant. Also 



C^{Hoi = {vi<S}:ie[l,k]} 



(3.25) 



gi{cti + d, ctk + rf) = cgi{ti, tk) + d 



(3.26) 
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for any constants c > and d. Then 

f{r]u...,rik,a) = ET{ ) (3.27) 

is nondecreasing in each rji when the other rjj and a are held constants. 



Proof of Theorem (4]. Due to (13.261) . we assume Yq has a mean and Yi has a mean 
Vi{= /^i "~/^o)- Let (f){x) be the pdf of A^(0, 1) and gu{y) be the pdf of a x^-distribution 
with u = I]f=o^j ~ + 1) degrees of freedom. Then 

fi7^,,...,7^,,a)= / ETi ^ ^ )^{x)g.iy)dxdy. 



For each fixed x and y, let 

Tx,y(Yi, Yfc) = T( 



which is non-decreasing in each % due to the monotonicity of T and gi^s. Therefore, 

the conditional distribution of T^^y for given x and y is stochastically nondecreasing 

in each //^(see Lemma 2 in Wu and Wang (2007)). Hence its conditional expectation 

Qi — x-^ — 5 Ofc — x-^ — 5 
ETx,y = ET{- ^ , # ) (3.28) 



is nondecreasing in each rji. So is /, the integral of (I3.28p . □ 

Remark 3. Each gi{Yi, Y^) = fii satisfies (13.261) and is nondecreasing in each Yj. 
We will use this to construct step-up tests in the next section. □ 

Remark 4. If define gi{Yi, ...,Yf:) = Yi for i G [^,k] and a sequence of 

for j G [1, A;], then gi and Tj^^ satisfy the conditions of Theorem SI Hsu and Berger's 
step-down tests (1999) claim A^, the MED, to be jo if Tf^^ = 1 but Tj^^^-^ = 0. □ 
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3.3 The construction of step-up tests 

We first construct a rejection region i?f with level a for each individual Hqi. 
Lemma 4 For a constant c, let 

i^le = { ^^"J°"^ >c}. (3.30) 

Then 

sup^^Ho^PfiiRic) = PnXRic), (3.31) 

where fi. = (/iq, /xi, /ifc) with /ii = ... = /ij = /io + 5 and /ij+i = ... = jJ-k = 
Therefore, for any a G (0, 1), Rl^^, with c = Ci^a, defines a level-a test for H^i, where 
Ci^a is the solution of 

PfiiRU = a. (3.32) 
Proof of Lemma [41. Let T = Ij^i . Then Lemma H] follows Theorem |H □ 



Remark 5. ci^a = ta,u^^/ni + l/rio due to (li = Yi when /x = /x^. □ 

Region satisfies property *), but not property **) in Section 2. To obtain 
more powerful multiple tests, we propose 

Theorem 5 For any integer i G [1, k] and for a sequence of nonnegative constants 
Ci through Ci, let 

i?.,...,,. = U}=i4,^ = > c,}. (3.33) 

Then 

SUp^^Ho,P^{Rcu...,c,) = P^i{Rc^,...,c,)- (3.34) 

Therefore, for any a G (0, 1), Rci,...,c,, with Ci = Ci^a given in Remark 5 and Ci 
determined iteratively by solving 
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P/xX^ci,...„cJ = «, (3.35) 
for i = 2, ...,k, defines a level-a test for i^oi- 

Proof of Theorem [5l Let T = Jo ^ . Then Theorem |5] foUows Theorem HI □ 

Theorem 6 Consider all hypotheses in C in l\3. 25\\ with the following testing proce- 
dure: 

assert Hai (or not H^i) if Ri =^ Rci,...,ci occurs (3.36) 

for any fixed a G (0, 1). Then the experimentwise error rate is at most a. i.e., the 
probability of making at least one incorrect assertion is at most a. 

Proof of Theorem [6l The proof is trivial if one notices that H^i is decreasing in 
i, and Ri is of level-a and is increasing in i. Then Theorem |6] follows the closed test 
procedure by Marcus, Peritz and Gabriel (1976). □ 

Remark 6. Region Ri is increasing in i. Then Ri satisfies properties *) and **). □ 

Remark 7. When the design is balanced, region Ri contains the set of {T^^ = 1}) 
on which Hsu and Berger's test (1999) claims the MED=1. Therefore, the proposed 
test is uniformly more powerful than Hsu and Berger's when the MED=1. □ 

Example l(continued). The sample standard deviation 5* = 7.751, to.o5,50 = 1.6 76. 
We compare the new step-up procedure with the step-up procedure SUIP, the step- 
down Williams procedure and step-down procedure SDIP in Tamhane et al (1996) 
and the DR method in Hsu and Berger (1999). For illustration, 6 = 6.5. From Table 1 
in Dunnett and Tamhane (1992) we have the critical values for the step-up procedure 
SUIP ci = 1.645, C2 = 1.933, C3 = 2.071, C4 = 2.165, C5 = 2.237, Cg = 2.294, C7 = 
2.342, cg = 2.382 (we treat df = 50 as df = oc). The SUIP procedure infers MED = 
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5. The Williams procedure has the t statistics: ti = — 1.810, )f:2 = —0.961,^3 = 
0.313 J4 = 1.899, ts = 5.788, te = 1 7 = ts = 9.334,^9 = 9.877. The Williams statistics 
are compared with the following critical values (taken from Williams (1971)) in a 
step-down manner: Ci = 1.675, C2 = 1.755, C3 = 1.780, C4 = 1.790, C5 = 1.795, Cg = 
1.800, C7 = 1.805, Cg = 1.805, cg = 1.810.. The Williams procedure infers AdED = 
4. By simulation with 1,0000 repetition, the critical values of the new statistic are 
ci = 0.968, C2 = 1.022, C3 = 1.046, C4 = 1.046, C5 = 1.034, ce = 1.043, C7 = 1.044, cg = 
1.047, Cg = 1.030, respectively. Also (/t^ - Yq - 6)/S for i = 1,2,3,4 given in fl333D 
are -1.045, -0.555, 0.181, and 1.097, respectively. Thus the new step-up procedure 
concludes MED = 4. So does Hsu and Berger (1999) 's DR method. 

4 Discussion. 

In this paper, we propose an alternative, SDMMSA, for the widely used PAVA. Al- 
though the two are equivalent, the SDMMSA is important by itself since it is easily 
coded and is well defined. Also a stochastic ordering of the estimators for the mono- 
tone normal means is established based on the SDMMSA. As one of its applications, 
a step-up test procedure is proposed to identify the MED. It strongly controls the 
experimentwise error rate, and is powerful to detect the MED, especially when the 
true MED is small. 
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TABLE 1. Sample dose-response data in Example 1 



Dosage Sample Yi SD Index /tj ij 



(mg/kg) 


size 




response 













6 


25.5 


2.6 





- 




- 


0.5 


6 


23.9 


4.0 


1 


23.9 


^7 


= 1 


1.0 


6 


27.7 


3.3 


2 


27.7 


^6 


= 2 


1.5 


6 


33.4 


2.3 


3 


33.4 


k 


= 3 


2.0 


6 


40.5 


10.5 


4 


40.5 




=-4 


2.5 


6 


57.9 


9.9 


5 


57.9 




= 5 


3.0 


6 


74.4 


14.6 


6 


73.77 




= 6 


3.5 


6 


73.4 


7.6 


7 


73.77 






4.0 


6 


73.5 


4.5 


8 


73.77 






4.5 


6 


76.2 


7.9 


9 


76.2 


k 


= 9 
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